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We perform a Bayesian model selection analysis in the the R-parity conserving MSSM to compare 
two different assumptions: whether the lightest neutralinos make all or only part of the cold dark 
matter. This corresponds to either imposing full WMAP relic density limits or just its upper 
bound for constraining the MSSM parameters. We consider several realisations of the MSSM, 
namely, three GUT-scale SUSY breaking scenarios with a handful of parameters corresponding to 
the CMSSM, anomaly mediation and the large volume string scenarios as well as the weak-scale 25- 
parameter phenomenological MSSM (pMSSM). The results give a data-based quantitative evidence 
for a multicomponent cold dark matter. The pMSSM posterior samples indicate that the choice of 
, imposing full WMAP limits or just its upper bound affects mostly the gaugino-higgsino content of 

■ the neutralino and, against naive expectations, essentially not any other sector. 

Astrophysical observations [l| indicate that more than 20% of the energy density of the universe is made of neutral 
and weakly interacting non-baryonic 'dark matter' (DM) (for reviews, see Refs. 041] )■ Neutrinos are known to 
make up a minor component of the total DM [5(. The bulk of the DM has to be cold in order to be consistent with 
structure formation in the universe. A cold DM (CDM) candidate should be massive and stable on cosmic time-scales. 
Oh' The neutralino lightest supersymmetric particle (LSP) in the minimal supersymmetric standard model (MSSM) with 
conserved R-parity is an excellent CDM candidate. Other CDM candidates from models beyond the standard model 
' include axions/axinos, gravitinos, Kaluza-Klein gravitons and stable moduli fields. 

In most of the MSSM literature the neutralino has been considered to be the only thermal relic contribution to 
the DM energy density. However this assumption may be too strong and unnecessary. In fact it is very difficult, 
i.e. requires much fine-tuning, to obtain the observed value of the CDM relic density with the above mentioned 
assumption. Typically the neutralino relic density either radically over-closes or under-closes the universe. For 
example bino-like neutralinos yield relic densities between 1 to 4 orders of magnitude more than the WMAP fit value, 
and it is known that in the minimal anomaly mediated symmetry breaking (mAMSB), which maps to a subset of the 
MSSM parameters space, neutralino co-annihilations are so important that the predicted relic densities are generically 
far below the cosmological and astrophysical fit value @ to the extent that an extra CDM component or some other 
modifications is unavoidable Q. Mixed axion-axino Q and/or stable moduli [9| could make up additional components 
to the neutralino CDM. In Rcfs. [Iol - fl2l ] multicomponent DM models were constructed to explain its indirect detection 
observations. 

The assumption that CDM are solely made of neutralinos and the alternative case where additional non-neutralino 
components are allowed are mutually exclusive hypotheses that could be compared against each other in light of the 
currently available indirect collider and astrophysical data. One could ask the question of which of the two alternative 
hypotheses has more support from current data relative to the other 71]. Answers to questions of this sort, which 
can be updated as new data become available, can be obtained via the Bayesian technique for model selection. 

Ba yes ian model selection techniques have been successfully applied in cosmology and astrophysics, see for instance 
Ref. [13| and references therein, and recently in particle physics phenomenology as well. For example the sign of 
the MSSM Higgs doublets mixing parameter, /i, is not yet fixed by any observation and is an important quantity in 
determining possible SUSY contribution to the muon anomalous magnetic moment, Sa^ [Til ]. In Ref. (lo1.ll6j] Bayesian 
technique were employed to select between the sign(/i) > and sign(fi) < branches of the constrained minimal 
supersymmetric model (CMSSM) and of the pMSSM. Another example which is conceptually different is that of 
Bayesian comparisons between different models as illustrated in Ref. [171 ] where various models for SUSY breaking 
mediation were compared with one another in light of current indirect collider and WMAP data. 

Bayesian approach to MSSM and DM phenomenology for the LHC have been used by various groups; see for 
instance Il8l-[2ll and references therein. In this article we compare the two CDM hypotheses in the context of the 
pMSSM [16|, [22|, (23[ and GUT-scale models for SUSY breaking mediations, using Bayes' theorem and current data, to 
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make statement about which of the hypotheses receive better support from the data. In the first hypothesis, 3~Coi the 
WMAP inferred CDM data is imposed on the pMSSM and other GUT-scale models parameters in a way compatible 
with the assumption that the whole CDM in the universe are solely made of neutralinos. For the second hypothesis, 
J£i, where a multicomponent CDM is assumed, pMSSM and other GUT-scale models points that give relic densities 
lower than WMAP central point are not penalised so that non-neutralino CDM components are allowed for. The way 
in which the WMAP CDM data is employed for the analysis is summarised in Fig. [lja) . 

The remaining part of the article is organised as follows. We first give a brief review of the Bayesian model selection 
technique followed by a description of how it is applied for the pMSSM and other GUT-scale models for the two CDM 
hypotheses. Then we compare and explain some differences between the sampled points, particularly in the gau gin g 
and neutralino sectors, from this analysis and the previous corresponding pMSSM analysis considered in Ref. |16| . 
Finally we explain the implications of the hypotheses for direct DM detection experiments before ending the article 
with summary and conclusions. 

Here we give definitions and set terminology that would allow for a self-contained discussion. Bayesian inference 
studies are centered around two quantities: model evidence and parameters posterior probability distributions. These 
are defined as follows. Consider a hypothesis (or model with all its assumptions) !H with parameters m. The a priori 
assumed form of values which the model parameters can take is encoded in the prior probability distribution of the 
model parameters given that the hypothesis is true, p(m\3~t). The posterior distribution of the parameters in light 
of the data considered and in the context of the hypothesis is represented by p(m\d, "K). The data, d, is a set of 
experimentally determined values for certain quantities (observables) which the model can predict. The support or 
evidence which the model obtains from the data is defined as the probability density of observing the data set given 
that the hypothesis is true and represented by Z = p{d\K). 

The evidence is calculated as 

Z = / p(d\m,3i)p(m\W) dm (1) 



where the integral is A^-dimensional, with N the dimension of the set of parameters m. p{d\m, J£) is called the 
likelihood, it is the probability of obtaining the data set d from certain model parameters m and is a function of x 2 
for the analyses presented in this article. Eq.([T]) is obtained directly from Bayes' theorem 

-r~\s art _ P(d\m,X)p{m\K) 

pWi) ■ (2) 

In the reign of searches for new physics [24j from upcoming LHC and other experimental data, Bayesian inference 
analysis would play an important role in selecting between various models based on the data. For such LHC^ 1 
problem and the Bayesian selection example presented in this article, evaluating the evidence is more important 
compared to computing parameters posterior probability distributions. In general one would not bother about finding 
a model posterior probability distribution if, to begin with, it is known that its evidence is poor. However, evaluating 
the evidence is computationally very expensive especially for multi-dimensional parameter spaces. In fact until the 
work in Ref. [25[ , the magnitude of the evidence were mostly disregarded and only the model posterior probability 
distributions were computed. For evaluating the evidence we use the MultiNest algorithm [26l ul\ which implements 
nested sampling technique [25j | and is very efficient in dealing with complex and multi-modal parameter spaces. 

In order to select between two hypotheses "Kq and "K-y one needs to compare their respective posterior probabilities 
given the data set d to be used for the comparison. The posterior probabilities are compared via Bayes' theorem as 

p{"Ki\d) = p(d|3fi)p(J{i) = Zi p(JCi) 

p(X \d) p(d\X Q )p(X ) Z p(X Q )' [ ' 

where p{'K\) / p{'Kq) is the ratio of prior probabilities for the two models usually set to unity with the assumption that 
the two hypothesis are a priori equally likely. Z\/Zq > 1 would mean that the data favours model J£i relative to 
J£o and vice versa if the ratio is less than one. Given various evidences to be compared beside one another, Jeffreys' 
scale [Hj], which gives a calibrated spectrum of significance for the relative strength between the evidences, 
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(4) 



is used. The Jeffreys' scale convention we employ is shown in Tab. fl] Values in the column labeled "Probability" are 
calculated as follows. Given that IKo and ?Ci are mutually exclusive and a priori equally likely then 

p(!K \d) +p(!Ki\d) = 1 andp(IKo) = j>(5Ci) (5) 
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|Alog e 


Z\ Probability 


Remark 


< 1.0 


< 0.750 


Inconclusive 


1.0 


0.750 


Weak Evidence 


2.5 


0.923 


Moderate Evidence 


5.0 


0.993 


Strong Evidence 



TABLE I: Jeffreys' scale for the interpretation of relative evidences. A log e Z > 5, e.g. 10, would indicates an irrefutable 
relative evidence. 



so that the posterior probabilities are given by 

p(^o \d) = — ^— and p(Xi \d) = — . (6) 

Next, we are going to apply the selection technique presented above to assess the mutually exclusive CDM hypotheses 
using current indirect collider and WMAP data. 

The pMSSM [l(| HH, US HH has unconstrained MSSM parameters defined at the weak-scale. In order to suppress 
sources of unobserved CP-violation and FCNC the parameters are set to be real with diagonal sfermion mass pa- 
rameters, rrif, and trilinear scalar couplings. In addition, first and second generation mass parameters are set to be 
degenerate. The soft SUSY breaking parameters that remain, out of the initially more than 100 parameters, make 
the pMSSM parameters set 

m= {tan/3, mL m 2 H ; M 12a ; m~ dgen , m 1 / 2nd9en ; At^.r.^e} (J) 

1 2 71,2,3,4,5 A, 2, 3,4,5 ' 

where the notations are as in Ref. [16[ and the integers 1,2,3,4,5 label third generation and degenerate first and 
second generation of the sfermion masses there. 

Cosmological and astrophysical fits to the standard cosmological constant plus CDM model imply a relic density 

m 

VcDAih 2 = 0.1143 ± 0.0034 (8) 

where h is the reduced Hubble constant [f| ■ Assuming that the thermal relic of the LSP account for the magnitude 
in Eq.©, the JCq hypothesised pMSSM predictions for the relic density are constrained to lie within the observed 
central value but with an inflated error to cater for uncertainties in theoretical predictions |30| 

VcDMh 2 =0.1143 ±0.02. (9) 

On the other hand, for "K\ only the upper half of the above constraint, Eq.©, is imposed such that predicted relic 
densities lower than the central value are not penalised as shown in Fig. 0Ja) . 

Apart from the CDM relic we also impose precision electroweak and B-physics indirect collider data for checking the 
neutralino CDM hypotheses. Various predictions, O, for the observables were obtained via S0FTSUSY2.0. 17 [31| for 
producing the MSSM spectrum; micrOMEGAs2 . 2 [32] for computing the neutralino CDM relic density, the branching 
ratio BR(B S fi + n~) and the anomalous magnetic moment of the muon (g — 2)^; SuperIso2 . }33j for predicting 
the Isospin asymmetry in the decays B — > K*j and BR(B — > s-f); and susyPOPE J3J, HH for computing precision 
observables that include the VF-boson mass my/, the effective leptonic mixing angle variable sin 2 @^ff, and the total 
Z-boson decay width, T z, at two loops in the dominant MSSM parameters. The experimentally determined central 
values, /i, and their corresponding uncertainties, a, of the observables 

O = {m w , sin 2 6 l e J f , T Z , Sa^, B$, A%, A 1 = A e , i?° c , A%, A b > c , (10) 
BR(B -> X s i), BR{B S -> n + fj,~), A _, Rbr{b u ^tv), Ram Bs , 

&>CDMh } 

make the data set (same as those in Ref. |l6|) for the analysis 

d={fi. l ,a i } (11) 

where i = 1, . . . , 19 label the individual observables. The central values and errors of the observables are summarised 
in Tab. [TTJ We assume that the observables are independent to form combined likelihoods 
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with full WMAP Gaussian constraint 

— with only WMAP upper bound constraint, flat prior pMSSM 
with only WMAP upper bound constraint, log prior pMSSM 
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FIG. 1: (a): The constraints put on pMSSM predicted values of SJcdmA 2 for the case where CDM are entirely due to the 
lightest neutralinos (solid line) or where a multicomponent CDM is allowed (dotted line) . (b) : The neutralino CDM relic density 
posterior distribution for Mo pMSSM hypothesis (solid line) and for Mi pMSSM hypothesis (dashed line.) For comparison the 
posterior from Ref. [l6l ] is also shown (dash-dotted line.) The posterior distributions show that generically the most preferred 
neutralino relic densities are centred around a point about two orders of magnitude less than the WMAP central value except 
when the single neutralino CDM scenario is imposed. This indicates that the single neutralino hypothesis is disfavoured. See 
Table II for a quantitave measure. 



Observable Constraint 



Th. Source Ex. Source 



Observable Constraint 



Th. Source Ex. Source 



mw [GeV] 


80.399 ± 0.027 


[34], [35] 


[36J 


T z [GeV] 


2.4952 ± 0.0025 


[37] 


[38] 


sin 2 e\j f 


0.2324 ± 0.0012 


[3JJ 


[38] 




(30.2 ± 9.0) x 10 10 


[39-42J 


[43-45] 


R l 


20.767 ± 0.025 


[37] 


[38] 


R b 


0.21629 ± 0.00066 


[3JJ 


[38] 


Rc 


0.1721 ±0.0030 




[38] 


A b FB 


0.0992 ± 0.0016 


[37] 


[38] 


A C FB 


0.0707 ± 0.035 


[37] 


[38] 



A = A e 
A b 
A c 

Br(B -+ X sl ) 
Br(B s fx + fl~] 
Ram Bs 

Rbt{B u ^tv) 

A - 

QcDAlh 2 



0.1513 ±0.0021 
0.923 ± 0.020 
0.670 ± 0.027 
(3.55 ± 0.42) x 10 4 
see Fig. [2] 
0.85 ±0.11 
1.26 ±0.41 
0.0375 ± 0.0289 
0.11 ±0.02 



[3D 


[38] 


[37] 


[38] 


[37] 


[38] 


[46-49J 


[50] 


[32, 51-53] 


[54] 


[55]~ 


[56] 


[57-59] 


[60-62 


[33] 


[63] 


[32, 51-53] 


[64] 



TABLE II: Summary for the central values and errors for the observables. 
quadrature to the experimental uncertainties quoted. 



Theoretical uncertainties have been added in 



Ml TT ™p[-(Oi-fi) 2 /2°f\ A ( v , „ ,TT ex P [-(0,-^) 2 /2q z 2 ] 
p(d|m,M ) =11 -=== - and p(d\m, JC X ) = I I ^ -=== 1 (12) 

i=i V 2 ^ i= i V 271 "^ 

where for "K\ the index j run over the different experimental observables (data) other than the CDM relic density; x 
represents the predicted value of neutralino CDM relic density; 



TM = ) 1 ^ y+ V 7 ™ 2 / 2 ), itx<y 

{) \exp[-(x-y) 2 /2 S 2 ]/( J/ ±^V 2 ), if x > y ' 1 J 

y = 0.1143 is the WMAP central value quoted in Eq.(|9]) and s = 0.02 the inflated error. L(x), shown in Fig. [TJa), 
is the likelihood function for Jto and %\ respectively corresponding to the regions where x < y and x > y. Next we 
present and discuss the results of the samplings in the paragraphs that follow. 

The Bayesian evidence values from the samplings are summarised in Tab. IIIII Let us first comment on the case 
of the pMSSM. For the setting of Ref. [23| , where only the linear prior was used, we obtain the following evidences: 
log e Z = 29.902±0.040, log e Z 1 = 33.120±0.028 and difference A log e Z = 3.218. Wider prior parameters ranges and 
more B-physics observables were added for the new samplings for this article and for those in Ref. [l6j relative to those 
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FIG. 2: The chi-squared contour used for the B-physics observable BR(B S — ¥ fi + /j,~). 



Models 


A log e Zfiat 


A log e Z log 


P{CDM = 


X?) 


pMSSM [23] 


3.2 




0.04 




pMSSM [16] 




2.7 


0.06 




CMSSM 


1.4 


2.5 


0.18 




mAMSB 


14.1 


14.8 


6.8 x HT 


7 


LVS 


3.6 


3.7 


0.02 





TABLE III: Alog e Zi values for the mutually exclusive CDM hypotheses in the pMSSM and some high scale models of SUSY 
breaking mediation. P(CDM = Xi) represents the probabilities that the CDM are solely made of neutralinos. For the details 
of the comparisons between the CMSSM, mAMSB and the LARGE volume scenario (LVS) [65[ models see Ref. [l7t ] . 

in Ref. [23j]; as such the absolute evidence values differ significantly. However, here only the difference in the evidence 
values are important for the purpose of this paper. A significant difference in A log e Z is not expected if for the flat 
prior the constraints _applied are as in Ref. (Tg in place of those in Ref. [23| from which the flat prior result quoted in 
Tab.|m]is obtained @. The log prior results in the pMSSM [3 are log e Z = 65.043 ±0.042, log e Zi = 67.761 ±0.030 
with A log e Z — 2.718. Translating the evidence ratios using Jeffreys' scale, Tab.|H implies that current indirect collider 
and cosmological data as employed here shows a significant evidence in support of multicomponent CDM hypothesis. 
The probability of an entirely neutralino-made CDM is around 0.04 or 0.06 respectively for the flat and log prior 
pMSSM. The pMSSM results imply an approximate prior independence since changing the log prior to a flat prior 
for the [l|| fit should not lead to a large difference in A log e Z. 

In the second part of Tab. IIII1 the difference in the evidence values for the two hypotheses Dio and JCi , is presented in 
the context of CMSSM, mAMSB and the LARGE volume scenarios (LVS) [HEEl]. Using the Jeffreys' scale Tab. Q] we 
can extract the following information: the CMSSM results are prior dependent ranging from being weak to moderate 
evidence for the multicomponent CDM hypothesis. The LVS is approximately prior independent with a clear-cut 
indication of the presence of significant evidence against the solely neutralino CDM hypothesis, Jfo- The mAMSB 
case is overwhelmingly prior independent and give an irrefutable evidence for multicomponent CDM hypothesis. The 
latter result can be considered as a check that proves that the Bayesian model selection technique is working since it 
is known [(J before hand that the mAMSB neutralino CDM cannot account for the WMAP fit, Eq.©, and hence a 
multicomponent CDM is necessary in that scenario. 

Note that from the evidence computations, the posterior samples of the models considered are obtained for free. 
For instance the posterior distributions of the predicted neutralino CDM for both hypotheses "Kq and 5f i are shown in 
Fig. [ljb). The distribution for J£i is centred at a value far less than the observed central value showing that addition 
non-neutralino CDM component (s) is (are) necessary in order to account for the WMAP data. Other analysis we 
performed from the posterior samples of the models are presented in the remaining parts of this article. We present 
and compare the difference between samples here for "K\ and from Ref. [161 ] for 3~Cq . We address mostly the neutralino 
and chargino sectors which get most affected by the choice of "Kq or IKi hypotheses. 

The nature of the neutralino determines the processes by which it annihilates and/or co- annihilates into standard 
model particles and therefore the main determining factor for its relic density. The neutralino mass is given by 
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\il> oT M N ^ + H.c. where tJj° T = (-ib, -iw 3 , H°, H$) and, 



Mk 



( Mi —mzcpsw mzspcw \ 

M 2 m z cpcw —mzspcw 

—mzcpsw mzcpcw — /j, 

\ m z sps w -mzspcw — A* J 



(14) 

c x = cosx and s x = sinx. The neutralino mass eigenstates are \i = ^ij^j where N is a unitary transformation that 
diagonalises M^. The neutralino mass eigenstate is a mixture of bino 6, wino w 3 and higgsinos Hi 2 

xi=N 11 b + N 12 w 3 +N 13 H°+N 14 Hl £ (A^) 2 = 1. (15) 

1=1,2,3,4 

The coefficients Nu with i = 1,2,3,4, and the neutralino masses are complicated functions of the soft terms 
Mi, M 2 , sign(ii), tan/3 etc (see e.g. [67j and references therein) however the following statements approximate 
some relations between the Nus and the parameters. When Mi -C min(M2, Nn dominates and the LSP is dom- 
inantly bino. Bino LSPs give relic densities too high, beyond the WMAP limits, in most regions of parameter space 
and therefore any pMSSM point in both Jfo and "Ki priors with such property would be excluded by the relic density 
constraint Eq.®. Next when M2 < min{M\, |/x|), N12 dominates so the the neutralino is dominantly wino and is 
quasi mass degenerate with the lightest chargino. This leads to strong chargino co-annihilations to the extent that 
the relic density is typically much smaller than the WMAP constraint. The pMSSM parameter points that have that 
type of neutralino would be excluded in the 'Kq prior scenario but not in For < min(Mi, M2), N13 and N14 
are of order one and the LSP is dominantly higgsinos-like. In this case the neutralino efficiently annihilates into top 
and weak gauge boson pairs. In addition the %?, X21 an d Xi are all approximately mass degenerate and higgsinos-like 
and hence there are more open channels for co-annihilations. Parameter points with this type of neutralino are also 
excluded in J£ prior. 

The preferred content-nature of the neutralinos, based on the data and CDM prior assumption considered, can be 
seen from the pMSSM posterior samples. We use a measure of the gauginos content Z g = lA^nl 2 + | iV22 1 2 , so 1 — Z g is 
approximately equal to unity or zero if the neutralino is mostly higgsino- or gaugino-like respectively. The posterior 
distributions of 1 — Z g for both J-Cq (solid line) and "K\ (dashed line) pMSSM priors are shown in Fig. [3Ja). It shows 
that the neutralino is mostly higgsino-like for the J£i prior pMSSM. The distribution in the J-Cq prior is bimodal 
with dominantly gaugino-like neutralinos and a sub-dominant mixed gaugino-higgsino peak. The drastic difference 
between the posteriors resulting from the two hypotheses is due to the fact that the CDM constraints imposed reject 
points that predict relic densities much less than the observed central value in the case of Dio . As result the posterior 
samples are dominantly gaugino-like but with sub-dominant mixed higgsino-gaugino neutralinos that allow moderate 
co-annihilations that keep the relic density within the Jfo allowed range. Jfi prior, on the other hand, allow and in fact 
prefer (as it turned out to be) channels that have very efficient co-annihilations which happen mostly for higgsino-like 
neutralinos. 

The neutralino relic density is inversely proportional to the thermally averaged magnitude of its annihilation and co- 
annihilation reaction cross-sections at early universe time. The content-nature of the neutralino determines which of 
the early universe reactions dominate over others. For instance, a certain neutralino composition lead to an acceptable 
magnitude for the thermally averaged cross-sections (which is proportional to inverse of the relic density) with 8% of the 
cross-sections coming from direct Xi annihilations into W- and Z-bosons, 63% from XrX® co-annihilations into quarks 
and 27% into leptons. Another example point has 46% and 48% of the thermally averaged cross-sections respectively 
coming from x^-m and x?-c co- annihilations. We keep record of this information for each of the pMSSM points visited 
during the nested sampling of the parameter space. The posterior distributions for the dominant annihilation and/or 
co-annihilation channels (or reactions) are shown in Fig. [3jb) for 3~C where the dominant channels are co-annihilations 
with sleptons and in Fig. [3Jc) for "Hi where co-annihilations with charginos and annihilations via chargino exchange 
dominate. 

The parameters posterior probability distributions are approximately the same for both Jfo and JCi pMSSM samples 
except in the neutralino and chargino sectors. This is mainly due to the difference in the properties of the neutralinos 
that lead to different co-annihilation processes which in turn control the consistency of the predicted relic densities 
with the imposed CDM constraints. The neutralino mass, chargino mass and the soft breaking parameters they depend 
on are shown in Fig. 0J The neutralino mass and tan j3 are not severely constrained by the relic density constraints 
and hence their posterior distributions do not drastically differ between the hypotheses. The other parameters and 
chargino masses, on the other hand, significantly differ. The reasons for this are explained below. There is also a 
wide difference in the posterior probability distributions of the SUSY breaking parameters Mi and M2 between the 
hypotheses since they determine the nature of the neutralino admixture which is tightly related to the relic density 
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FIG. 3: (a) Neutralino's gaugino-higgsino admixtures for "Kq (solid line) and for K 1 (dashed line) log prior pMSSM CDM 
hypotheses. The neutralino is almost always purely higgsino (peak around zero) for the Jii hypothesis due to a single efficient 
dominant co-annihilation that leads to posteriors with much lower (than WMAP central value) relic densities. For the "Ko 
hypothesis, however, there are mainly two dominant (co-)annihilation channels that keep the relic densities within the imposed 
WMAP limit. With Z g — |iVu| 2 + | TVaa | 2 , peaks to the right-hand side of the abscissa indicate higgsino domination, (b) !Ko 
prior pMSSM posterior sample distribution for the early universe neutralino annihilation (solid line) and co-annihilations whose 
cross-sections dominate in determining the present neutralino relic density. Neutralino-sleptons co-annihilation is dominant in 
this scenario since the neutralinos are mostly gaugino-like as shown in plot (a) above, (c) Same as in (b) above but for !Hi 
prior pMSSM posterior samples. Here chargino-neutralino co-annihilations dominate. 



predictions via the preferred co-annihilation cross-sections. The /i parameter on the other hand is not as widely 
different since it is mostly controlled by the requirement of radiative electroweak symmetry breaking (EWSB) than 
by the CDM constraint. 

The nature of the Higgs doublets mixing parameter, /i, and mass parameters rriH 1 2 are controlled at leading order 
by the tree level requirement that 



1 2 rr? H — m 2 H2 tan 2 /3 
2 Wz= 'tan 2 /?-! 



m|= a \ ' — fJ- (16) 
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2 m, z _ —m H2 — /i with moderate to high tan /? values. That is for a fixed Z-boson 
mass, TOjy 2 and /x 2 would be roughly the same since renormalisation group (RG) running of the soft SUSY breaking 
parameters take m^r 2 to very large (relative to the gaugino mass parameters) negative values at the EWSB scale. 
This explains the shape of the diagrams for /j, 2 and vr? H parameters in Fig. SJ The 5£i pMSSM fj, 2 posterior samples 
(shown in dashed line in the figure) are less than those in the "Kq prior because in the former case the neutralinos are 
mostly higgsino-like compared to the dominantly gaugino-like neutralinos in the latter. The details of the neutralino 
gaugino-higgsino content or admixture is shown in Fig. [3ja). 

The posterior distributions for the gaugino mass parameters Mi and M2 widely differs between the two CDM prior 
assumptions since these parameters control the nature of the neutralino gaugino-higgsino admixtures. We are going 
to address the Mi and M2 parameters structure for each of the priors, starting with %i. The scenario in the J£i prior- 
is similar to the so-called focus point region in CMSSM where the m# 2 RG running is independent of the gaugino 
and trilinear parameters, see e.g. (68j . As a result Mi and M2 would in principle remain unconstrained by the 
EWSB requirement. This indeed remains the case for the Mi parameter as shown in the corresponding plot in Fig. 2] 
However M2 is constrained from the chargino sector since neutralino-chargino co-annihilations requires fn % ± ~ m x1 
which in turn imposes Mf ~ /1 2 . The M2 posterior distribution is similar to that of the /x parameter as shown in 



8 





II 




1 


I 





-4-2 O 2 4 





0.5 1 1 .5 2 

m o 






\\ 


if 




1 









20 40 60 

tan f3 




FIG. 4: Marginalised one-dimensional posterior distributions for selected pMSSM parameters that are most affected by the 
choice between the two hypotheses "Kq and JCi. Except for tan/3, all numbers in the horizontal axis are in TeV units. All the 
dashed lines are for the prior hypothesis Jfi where neutralinos form only part of multicomponent CDM. Solid lines are for the 
prior hypothesis IKq where neutralinos make up all of the CDM. 



Fig. S] (dashed lines). Now turning to the 5£o prior: again moderate tan/3 implies m 2 H2 ~ /Lt 2 but here the neutralinos 
are mostly gauginos-like with Mf > y 2 > M 2 and therefore different from the "K\ focus point-like feature. Mi is 
constrained to be small by co-annihilation scenarios requiring mass of the neutralino, controlled by M\ < M 2 , to be 
approximately degenerate with the sleptons it is co-annihilating with. In summary M 2 controls the mass of the mostly 
higgsino neutralinos via neutralino-chargino co-annihilations in JCi prior where the neutralinos make only part of the 
CDM energy of the universe. On the other hand if neutralinos were to make all the CDM then M\ < M 2 controls its 
mass via co-annihilation with with sleptons. 

We next turn to the chargino masses' posterior distributions. Co-annihilations with charginos constrain the chargino 
masses to be approximately equals those of the neutralinos. This is more so in the case of the Jf 1 prior (dashed lines 
in Fig. 2]) than in "Kq prior (solid lines) pMSSM since in the latter co-annihilations with sleptons dominate and hence 
the chargino masses can be much larger than the neutralino masses. The relation between m ± and m ± for each of 

Xi X2 

the priors can be explained as follows. The masses are controlled by the M 2 and fi parameters 

1 " 



xf,xf 2 



|M 2 | 2 + H 2 + 2m 2 w T J(\M 2 \ 2 + H 2 + 2m 2 w ) 2 - %M 2 - m 2 w sin2/3| 



(17) 



When Mf ~ /i 2 , as is the case for the 3ii prior pMSSM, then the following approximation holds true (|M 2 
iwf ~ %M 2 - m 2 w 



fi\ 2 + 2m 2 v ) 2 ~ 4|/iA/ 2 — m 2 v sin2/3| 2 . This in turn would imply, from Eq. dTTl) . that m 2 ± ~ m 2 ± as can be seen in 



X 2 Xi 

the corresponding plots in Fig. 2J In a similar manner when Mf >> /i 2 as is the case for !Ko prior pMSSM, then 
\/(|M 2 | 2 + \fi\ 2 + 2m 2 v ) 2 - 4|/uM 2 - m 2 w sin2/3| 2 - M 2 2 . This way m 2 ± - \y? and m 2 ± ~ M| explain the structure 

Xi X2 

of the corresponding "Kq distributions in Fig. |U 

Recapitulating, we have analysed the two mutually exclusive CDM hypotheses; whether it is solely made of neutrali- 
nos, Jfo, or whether a multicomponent scenario must be allowed, JCi. We applied the Bayesian approach to inference 
and used the current indirect collider and cosmological data to select between these hypotheses within the context 
of the pMSSM where 20 phenomenologically viable MSSM parameters are simultaneously varied at the electroweak 
energy scale. The pMSSM setting is unbiased and free of theoretical assumptions and uncertainties from sources 
of SUSY breaking, mediation mechanisms and RG running. We also applied the Bayesian technique to GUT-scale 
defined models of SUSY mediation mechanisms. The results shown in Tab. IIIII give a data-based evidence that a single 
neutralino CDM hypothesis is unrealistic and improbable. Based on the set of data employed for the analysis, the 
probability that neutralino make all the CDM for the pMSSM, mAMSB, CMSSM, and LVS are about 0.06 at most, 
6.8 x 10~ 7 , 0.18, and 0.02 respectively. The results in Tab. IIIII and Jeffrey's scale (shown in Tab. [TJ) for interpreting 
the evidence values shows that prior dependence of the results are strongest in the CMSSM, despite the fact that 
it has fewer number of parameters compared to the pMSSM. CMSSM has a weak-to-moderate evidence against J£o, 
the solely neutralino made CDM hypothesis. The pMSSM and LVS results have moderate-to-sub-strong evidence 
against IK - LVS result is approximately prior independent. The mAMSB model have an exceptionally strong and 
prior independent evidence against J£q- Based on these we therefore conclude that MSSM phenomenology studies 
should not neglect the mixed CDM possibility. 
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Let us, here, re-emphasise that prior dependence is a positive feature of Bayesian methods because it shows when 
the data and/or the physics under consideration are strong or well-defined enough to be free of ambiguous predictions 
and discriminations between alternative hypotheses. Prior (in)dependence is not a feature simply determined by the 
number of parameters in the sense that models with more parameters necessarily exhibit more prior dependence in 
their predictions. We hope that the results presented in this paper give good enough illustration of the mentioned 
observations. First, consider the mAMSB results in Tab. lIIII The results are prior independent because the physics is 
very well-pinned to the extent that even the not-that-strong indirect collider and cosmological data, Eq. (|lll) , employed 
is good enough to give such exceptionally unambiguous results. Prior dependence can be completely absent if the 
data employed is strong enough. Secondly, the presence of relatively stronger prior dependence in the model with 
less number of parameters, CMSSM, compared to that with more parameters, the pMSSM, illustrates that prior 
dependence is not simply a feature directly derived from the number of model parameters but is only a feature 
indicating the strength of the data employed and/or the physics under consideration. 

It is worth mentioning that it may seem that the analysis is not fair to the !Kq prior pMSSM, where neutralinos 
are assumed to form all CDM. Not fair in the sense that any additional CDM candidate would be accompanied by its 
own parameters describing its mass and couplings. It would seem that these additional parameters would dilute the 
prior and in turn reduce the evidence for the scenario. We emphasise that this is necessarily not going to be the case 
because the additional DM parameters would make the model under consideration something completely different 
from the pMSSM. As such the comparison would then be between a pMSSM and a non-pMSSM. And that is not 
what we considered in this article. The comparison carried out is strictly between same pMSSMs with two distinct 
and mutually exclusive hypotheses: in the first, ?£q, the CDM is assumed to be solely made of the neutralino while 
in the second, Jti , the neutralino forms only part of the CDM with the remaining to be accounted for by some other 
physics external and perpendicular to the pMSSM. 

The posterior samples resulting from the analysis with "Kq, the assumption that the neutralinos make all CDM, and 
JCi, where the neutralinos form only part of the CDM, hypotheses differ drastically in the neutralino gaugino-higgsino 
content and neutralino/gauginos sector parameters. The dark matter relic density constraint in 3~Cq dis-favours pMSSM 
points with higgsino-like or mixed gaugino-higgsino neutralinos that have efficient co-annihilation channels. However 
relaxing the constraint to Jf i disfavoured pMSSM points become the preferred ones. Results in both this and previous 
pMSSM work indicate that, in light of the indirect collider and cosmological data considered, neutralino-chargino co- 
annihilations dominate at the early universe time before the neutralinos decouple from thermal equilibrium. It would 
be interesting to see what implications this feature would have for MSSM phenomenology at colliders. 

Our results provide non-trivial quantitative evidence about the composition of DM in the MSSM, even after having 
introduced large number of parameters as in the pMSSM. We have assumed thermal dark matter and R-parity 
conserving MSSM. Modification of these assumptions or further extensions of the MSSM may modify our results. 
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the University of Cambridge for direct access to the super- computers: COSMOS from the Department of Applied 
Mathematics and Theoretical Physics (DAMTP) and the Darwin cluster from the High Performance Computing 
Service (HPCS). 



[1] E. Komatsu et al. (2010), 1001.4538. 

[2] G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rept. 267, 195 (1996), hep-ph/9506380. 
[3] G. Bertone, D. Hooper, and J. Silk, Phys. Rept. 405, 279 (2005), hep-ph/0404175. 
[4] K. Garrett and G. Duda (2010), 1006.2483. 

[5] E. Komatsu et al. (WMAP), Astrophys. J. Suppl. 180, 330 (2009), 0803.0547. 

[6] G. F. Giudice, M. A. Luty, H. Murayama, and R. Rattazzi, JHEP 12, 027 (1998), hep-ph/9810442. 

[7] H. Baer, R. Dermisek, S. Rajagopalan, and H. Summy (2010), 1004.3297. 

[8] H. Baer, A. D. Box, and H. Summy (2010), 1005.2215. 

[9] J. P. Conlon and F. Quevedo, JCAP 0708, 019 (2007), 0705.3460. 
[10] K. M. Zurek, Phys. Rev. D79, 115002 (2009), 0811.4429. 
[11] I. Cholis and N. Weiner (2009), 0911.4954. 
[12] D. Feldman, Z. Liu, P. Nath, and G. Peim (2010), 1004.0649. 
[13] R. Trotta, Mon. Not. Roy. Astron. Soc. 378, 72 (2007), astro-ph/0504022. 
[14] P. von Weitershausen, M. Schafer, H. Stockinger-Kim, and D. Stockinger (2010), 1003.5820. 
[15] F. Feroz et al., JHEP 10, 064 (2008), 0807.4512. 

[16] S. S. AbdusSalam, B. C. Allanach, F. Quevedo, F. Feroz, and M. Hobson, Phys. Rev. D81, 095012 (2010), 0904.2548. 
[17] S. S. AbdusSalam, B. C. Allanach, M. J. Dolan, F. Feroz, and M. P. Hobson, Phys. Rev. D80, 035017 (2009), 0906.0957. 



10 



[18] B. C. Allanach, Phys. Lett. B635, 123 (2006), hep-ph/0601089. 

[19] M. E. Cabrera, J. A. Casas, and R. Ruiz de Austri, JHEP 03, 075 (2009), 0812.0536. 

[20] A. R. Raklev and M. J. White (2009), 0911.1986. 

[21] M. J. White and F. Feroz, JHEP 07, 064 (2010), 1002.1922. 

[22] A. Djouadi et al. (MSSM Working Group) (1998), hep-ph/9901246. 

[23] S. S. AbdusSalam, AIP Conf. Proc. 1078, 297 (2009), 0809.0284. 

[24] P. Nath et al., Nucl. Phys. Proc. Suppl. 200-202, 185 (2010), 1001.2693. 

[25] J. Skilling, in American Institute of Physics Conference Series, edited by R. Fischer, R. Preuss, and U. V. Toussaint (2004), 

pp. 395-405, URL http: //www . inference .phy . cam. ac .uk/bayesys/ 
[26] F. Feroz, M. P. Hobson, and M. Bridges, Monthly Notices of the Royal Astronomical Society 398, 4, 1601 (2008), 0809.3437. 
[27] F. Feroz and M. P. Hobson, Monthly Notices of the Royal Astronomical Society 384, 2, 449 (2007), 0704.3704. 
[28] H. Jeffreys, Theory of probability, 3rd edn (Oxford University Press, 1961). 
[29] C. F. Berger, J. S. Gainer, J. L. Hewett, and T. G. Rizzo, JHEP 02, 023 (2009), 0812.0980. 
[30] N. Baro, F. Boudjema, and A. Semenov, Phys. Lett. B660, 550 (2008), 0710.1821. 
[31] B. C. Allanach, Comput. Phys. Commun. 143, 305 (2002), hep-ph/0104145. 
[32] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov (2008), 0803.2360. 
[33] F. Mahmoudi, Comput. Phys. Commun. 178, 745 (2008), 0710.2067. 

[34] S. Heinemeyer, W. Hollik, D. Stockinger, A. M. Weber, and G. Weiglein, JHEP 08, 052 (2006), hep-ph/0604147. 
[35] S. Heinemeyer, W. Hollik, A. M. Weber, and G. Weiglein, JHEP 04, 039 (2008), 0710.2972. 
[36] M. Verzocchi, in talk at ICHEP 2008 (2008, Philadelphia, USA). 

[37] A. M. W. et al., SUSY-POPE (Precision Observables Precisely Evaluated) (in preparation, ????). 

[38] Phys. Rept. 427, 257 (2006), hep-ex/0509008. 

[39] T. Moroi, Phys. Rev. D53, 6565 (1996), hep-ph/9512396. 

[40] G. Degrassi and G. F. Giudice, Phys. Rev. D58, 053007 (1998), hep-ph/9803384. 

[41] S. Heinemeyer, D. Stockinger, and G. Weiglein, Nucl. Phys. B690, 62 (2004), hep-ph/0312264. 

[42] S. Heinemeyer, D. Stockinger, and G. Weiglein, Nucl. Phys. B699, 103 (2004), hep-ph/0405255. 

[43] G. W. Bennett et al. (Muon G-2), Phys. Rev. D73, 072003 (2006), hep-ex/0602035. 

[44] M. Davier, Nucl. Phys. Proc. Suppl. 169, 288 (2007), hep-ph/0701163. 

[45] D. W. Hertzog, J. P. Miller, E. de Rafael, B. Lee Roberts, and D. Stockinger (2007), 0705.4617. 
[46] M. Misiak et al., Phys. Rev. Lett. 98, 022002 (2007), hep-ph/0609232. 
[47] M. Misiak (2006), hep-ph/0609289. 

[48] M. Misiak and M. Steinhauser, Nucl. Phys. B764, 62 (2007), hep-ph/0609241. 
[49] T. Becher and M. Neubert, Phys. Rev. Lett. 98, 022003 (2007), hep-ph/0610067. 
[50] E. Barberio et al. (Heavy Flavor Averaging Group (HFAG)) (2007), 0704.3575. 

[51] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 176, 367 (2007), hep-ph/0607059. 

[52] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 174, 577 (2006), hep-ph/0405253. 

[53] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 149, 103 (2002), hep-ph/0112278. 

[54] T. Aaltonen et al. (CDF), Phys. Rev. Lett. 100, 101802 (2008), 0712.1708. 

[55] M. Bona et al. (UTfit), JHEP 10, 081 (2006), hep-ph/0606167. 

[56] A. Abulencia et al. (CDF), Phys. Rev. Lett. 97, 242003 (2006), hep-ex/0609040. 

[57] G. fsidori and P. Paradisi, Phys. Lett. B639, 499 (2006), hep-ph/0605012. 

[58] G. fsidori, F. Mescia, P. Paradisi, and D. Temes, Phys. Rev. D75, 115019 (2007), hep-ph/0703035. 

[59] A. G. Akeroyd and S. Recksiegel, J. Phys. G29, 2311 (2003), hep-ph/0306037. 

[60] B. Aubert et al. (BABAR), Phys. Rev. Lett. 95, 041804 (2005), hep-ex/0407038. 

[61] P. Chang, in talk at ICHEP 2008 (2008, Philadelphia, USA). 

[62] A. Gray et al. (HPQCD), Phys. Rev. Lett. 95, 212001 (2005), hep-lat/0507015. 

[63] C. Amsler et al. (Particle Data Group), Phys. Lett. B667, 1 (2008). 

[64] E. Komatsu et al. (WMAP), Astrophys. J. Suppl. 180, 330 (2009), 0803.0547. 

[65] J. P. Conlon, S. S. Abdussalam, F. Quevedo, and K. Suruliz, JHEP 01, 032 (2007), hep-th/0610129. 
[66] S. S. AbdusSalam, J. P. Conlon, F. Quevedo, and K. Suruliz, JHEP 12, 036 (2007), 0709.0221. 
[67] M. M. El Kheishen, A. A. Aboshousha, and A. A. Shank, Phys. Rev. D45, 4345 (1992). 
[68] J. L. Feng, K. T. Matchev, and T. Moroi, Phys. Rev. D61, 075005 (2000), hep-ph/9909334. 
[69] J. D. Lewin and P. F. Smith, Astroparticle Physics 6, 87 (1996), ISSN 0927-6505. 
[70] O. Buchmueller et al., JHEP 09, 117 (2008), 0808.4128. 

[71] This question is different from attempts as in Ref. seeking to estimate the effects, on preferred CMSSM regions, of 

changing the error bars on the (solely neutralino) CDM relic density constraint from WMAP. 
[72] Here we use the five years WMAP values in order to properly compare with Ref. [16j. Updating this to the seven years 

WMAP value would not change the main results and conclusions of this work. 
[73] The main difference between the pMSSM analysis in Ref. [2^] and Ref.[l(| is that the latter have additional set of electroweak 

observables from the code susyPOPE [34l . [35| and have extended parameter ranges up to 4 TeV compared to 2 TeV for the 

former. 



